Three-dimensional matter-wave vortices in optical lattices 
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We predict the existence of spatially localized nontrivial vortex states of a Bose-Einstein con- 
densate with repulsive atomic interaction confined by a three-dimensional optical lattice. Such 
vortex-like structures include planar vortices, their co- and counter-rotating bound states, and dis- 
tinctly three-dimensional non-planar vortex states. We demonstrate numerically that many of these 
vortex structures are remarkably robust, and therefore can be generated and observed in experiment. 
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The importance of three-dimensional (3D) vortex con- 
figurations in physics was recognized long ago by Sir 
W. Thompson when he suggested that robust vortex 
rings can be employed to understand the structure of 
the atom 0] . Since then topologically stable vortex lines 
have been investigated in many physical contexts, includ- 
ing cosmology, hydrodynamics, optics, and condensed 
matter physics. New experimental opportunities for the 
study of vortices have appeared in the rapidly developing 
field of Bose-Einstein condensates (BECs) |2( as conden- 
sates provide an excellent test system for the study of 
nucleation, dynamics, and the interaction of vortex lines. 

One of the outstanding physical problems is the behav- 
ior of vortices in systems with reduced symmetry, such 
as periodic structures. BECs loaded into periodic po- 
tentials of optical lattices provide an excellent example 
of such a system, being relatively easy to manipulate ex- 
perimentally. In a two-dimensional (2D) geometry, stable 
vortices localized in an optical lattice were predicted to 
exist in the form of matter-wave gap vortices 0- These 
vortices are spatially localized states of the condensate 
with a phase circulation that exist due to the interplay 
of the repulsive nonlinearity and the anomalous diffrac- 
tion of the matter waves due to their Bragg scattering in 
the lattice. Although 2D (planar) vortices have also been 
investigated in 3D discrete lattice models |^|, nothing is 
known about nonlinear dynamics of the vortex states in 
a fully 3D continuous lattice geometry. 

In this Letter we not only show that planar 2D matter- 
wave gap vortices persist in the presence of the third 
lattice dimension, but also predict the existence of com- 
pletely new vortex states localized inside the lattice. In 
particular we show that the circulation of particles can 
break the traditional 2D vortex symmetry and form 
states which circulate simultaneously in multiple sym- 
metry planes. Furthermore we show that the 3D lattice 
provides a unique opportunity to study the coherent in- 
teraction of multiple vortices, and demonstrate that an 
interacting stack of vortices containing an extended vor- 
tex line, which is known to be unstable in homogeneous 
systems, is stabilized by the lattice. As 3D optical lat- 
tices provide improved control over condensate excita- 



tions 5], there is promise for the future experimental 
study of these highly nontrivial vortex phenomena. 

We model the macroscopic dynamics of the Bose- 
Einstein condensate in a 3D optical lattice by the rescaled 
mean-field Gross-Pitaevskii (GP) equation: 



•<9* A T 
i— + A* 

at 



T/(r)^ + 7|*|* = 0, 



(1) 



where ^> is the complex matter-wave wavefunction, r = 
(x, y, z), A = d 2 / dx 2 + d 2 / dy 2 + d 2 / dz 2 , and 7 = -1 for 
repulsive atomic interactions. The characteristic spatial 
scale is — d/n, where d = A/2 is the lattice period, 
and A is the wavelength of the light producing the lattice. 
We consider a cubic lattice potential: 



V(r) = V Q (sin 2 



sin 2 y - 



sin 2 z) 



(2) 



with the depth Vq scaled by the lattice recoil energy 
El = h 2 /2ma 2 L . The characteristic time scale is given 
by u>l = hjEjj- We have neglected any additional trap- 
ping potential as we are interested in spatial localization 
due to the optical lattice and nonlinearity. 

We look for the stationary solutions of Eq. @J using the 
ansatz ^(r, t) = ip(r) exp(-int), where /i is the chemical 
potential. At low atomic densities, the condensate dy- 
namics in the lattice is described by the linearized Eq. (JTJ 
where the last term is neglected. Due to the periodicity of 
the potential, the corresponding stationary states can be 
expressed as a superposition of the Floquet-Bloch modes 
satisfying the condition, ipk(r) = i?k(r)exp(ikr), where 
^k(r) = -Bk(r + d), and k is the quasi- momentum within 
the first Brillouin zone of the lattice (see Fig. QJ. Mat- 
ter waves can propagate through the lattice if k is real- 
valued, i.e. when the chemical potential of the Floquet- 
Bloch linear states belongs to one of the bands (shaded 
in Fig. nj. In the 3D potential (JSJ), there always exists 
a semi-infinite spectral gap extending to /i — > —00, and 
additional gaps appear only when the potential depth Vq 
exceeds the threshold value Vth ~ 2.24. In what follows 
we choose the value Vq — 6. 

At higher atomic densities, the effect of repulsive atom- 
atom interactions described by the nonlinear term in 
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FIG. 1: Matter-wave spectrum ^i(k) in a three-dimensional 
optical lattice @ at Vb = 6. Shaded and open areas show 
bands and gaps, respectively. Inset shows the first Brillouin 
zone and the high-symmetry points in the reciprocal lattice. 

Eq. Q becomes important. In free space, repulsive non- 
linearity tends to accelerate the spreading of the conden- 
sate. However, in a lattice the effective mass of the con- 
densate becomes negative near the lower gap edges 
allowing for the formation of localized states inside the 
spectral gaps |6|. These nonlinear states - gap solitons 
- are localized in all three dimensions. Figure [3 shows 
an example of a gap soliton, and the dependence of its 
norm, N = J \ ^>\ 2 dr, on the chemical potential within 
the first gap calculated using a numerical functional min- 
imization technique 7] . Despite a nontrivial phase of the 
wavefunction [see Fig.^c)], characteristic of linear Bloch 
wave tails, there is no particle flow in the soliton. 

Near the lower gap edge the soliton spreads over many 
lattice sites and may be described by the free-space GP 
equation for an effective envelope of the correspond- 
ing Bloch state Q. In this limit, the 3D gap soliton 
is expected to be unstable according to the Vakhitov- 
Kolokolov criterion since dN/dfi < [see Fig. Ufa)]. 
Such an instability may result in the transformation of a 
weakly localized atomic cloud into a highly confined state 
deeper inside the gap [such as shown in Fig.|2tb,c)] , where 
dN/dfi > 0. We have investigated the dynamical stabil- 
ity of the strongly localized gap solitons numerically by 
introducing 5% random phase and amplitude perturba- 
tions. No evidence of instability can be seen over the 
evolution time t = 10 3 , which for 87 Rb in an optical lat- 
tice created with a 0.25/im period is equivalent to 17ms. 

Spatially localized stationary states of the repulsive 
condensate with a phase circulation can also exist in the 
gaps of the linear spectrum. These are localized vortices 
with stationary density distribution, but nonvanishing 
particle flow around a closed path, which can be charac- 
terized by the local current density, j = — 2Im(>I'V r , I'*). 
The gap vortex states in a lattice can be thought of as 
clusters of gap solitons ip s having different positions on 
the lattice r m and phases (f> m . The vortex flow between 
such solitons occurs due to a nontrivial phase structure 
of the full matter-wave function which is approximately 




FIG. 2: (color online) Gap soliton in a 3D optical lattice with 
Vb = 6. (a) Norm of the wavefunction N vs. ri inside the gap; 
point C (at /i = 7.9) corresponds to the state shown in (b) 
and (c). Line at left band edge calculated in free-space GP 
limit. Dotted line an interpolation, (b) An isosurface of the 
soliton density, |V>| 2 , at the level l^l 2 = 0.2, with phase in the 
z = plane (semi-transparent) shown below, (c) The cuts of 
ij) (dashed) and \4>\ 2 (solid) along the line y = z = 0. Thin 
dashed line corresponds to |t/>| 2 = 0.2. Shaded areas show the 
lattice minima (parts of the potential less than Vb/2). 

found as ^ ~ X)£f =1 ^ a (r - r m )exp(z0 m - ifit). For- 
mation of a stationary vortex is possible when all the 
incoming and outgoing flows are balanced and hence [tj : 

M 

^ c nm sm(c/) m -<j) n ) = 0. (3) 

m—l 

Each term in this sum defines the flow between gap soli- 
tons numbered n and m, which depends on the rela- 
tive phases and the nonlinear coupling constants c nm = 
c mn = J J ipU r ~ r n)ips{r - r m )dr/ f ^(r)dr. 

The sign of the coupling depends on the sign of the 
overlap between the wavefunctions of the gap solitons and 
the sign of the nonlincarity. The former may be positive 
or negative due to the oscillating tails of the gap solitons 
[Fig. Hfc)]. As the latter is negative (7 = —1), positive 
coupling occurs when the overlapping parts of the wave- 
functions are 7T out-of-phase. Solitons in the cluster can 
be positioned in the "nearest neighbor" configuration, 
whereby they are separated by a single lattice period in 
the x, y or z direction, or in a "next neighbor" configura- 
tion, whereby the solitons are placed in-plane, diagonally 
across a lattice cell (e.g. in the direction r — ► M in Fig. 
1). The "nearest neighbor" coupling is positive, whereas 
the "next neighbor" coupling is negative. Importantly, 
for a positive coupling, the flow is directed from a soliton 
with phase cf> m = to one with tt/2 to ir etc., while for 
negative coupling the flow is reversed. Below we discuss 
the physical origins of this reversal. 
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FIG. 3: (color online) Examples of planar (a) on-site and (b) 
off-site vortices (Vb = 6, fx = 7.9). The orientation in (a) is 
the (x — y) plane, while in (b) it is the (y — z) plane. Arrows 
show directions of particle flow. 



A 3D lattice can support vortex states of different ge- 
ometries. First and the most obvious type are the gap 
vortices localized in a plane which can be considered as 
a generalization of 2D gap vortices Such 2D states 
were shown to have two basic symmetries relative to the 
lattice sites: on-site and off-site vortices yj. The off-site 
vortex has a lattice maximum at its center and, therefore, 
it may consist of four lobes (atomic density maxima) in 
a nearest neighbor configuration [Fig. Eta)]. As the near- 
est neighbor coupling is positive, the flow circulation is 
defined in the direction — > 7r/2 — > tt — > 37r/2 — > of 
the phases of the composite solitons. The on-site vor- 
tices are centered on a lattice site and consist of four 
lobes in the next neighbor configuration [Fig.Etb)]. With 
the nearest neighbor coupling normalized to +1, the on- 
site vortex has an inter-lobe next neighbor coupling of 
—0.11, and the particle flow is opposite to that expected 
from the composite soliton phases. At the physical level, 
as particle flow is always directed along the phase gradi- 
ent this seemingly anomalous flow actually arises through 
the interference of multiple phase singularities producing 
a counter-flow opposite to that of the central singularity. 

We have examined the dynamical stability of the pla- 
nar vortices shown in Figs.EI a jb) numerically and found 
that they are stable, at least up to t = 10 3 . An im- 
portant consequence of the 3D geometry is that we may 
orient these planar vortices in any of the three symme- 
try planes of the lattice. Additional planes exist in the 
three-dimensional lattice, however they are asymmetric 
with respect to the inter-site coupling, and we expect 
such planes to host asymmetric planar vortices . 

A new feature of 3D optical lattices is the possibil- 
ity of studying the coherent interaction of multiple pla- 
nar vortices. From the analysis of the discrete equations 
©, we have identified three basic types of coupled pairs 
of planar vortices: in-phase, out-of-phase, and hybrid- 
phase. These may be described by the phases of the 
vortex lobes. The in-phase pair has a phase structure, in 
multiples of n, corresponding to ([0,0], [5, \], [1, 1], [§, §]) 
where we have used the notation of [vortexi, vortex2] for 
each of the four lobes. The out-of-phase pair has a phase 




FIG. 4: (color online) Examples of interacting vortex soli- 
tons at Vo = 6, /i = 7.9. (a) In-phase co-rotating vortex pair, 
(b) Counter-rotating vortex pair. Lines correspond to pri- 
mary vortex lines, with different shadings showing opposite 
topological charges. 

of ([0, 1], [5, §], [1, 0], [f , \}), while the hybrid vortex pair 
has a phase structure ([0, 0], [|, |], [1, 1], [I, 5]). In an 
equivalent hybrid state the second vortex is shifted in 
phase by n. The particle flows of the in-phase and out-of- 
phase vortices are co-rotating, and hence the pair shares 
a singular vortex line at the center [Fig. Ufa)]. In con- 
trast, the hybrid vortex states are equivalent to interact- 
ing counter-rotating vortices. Their vortex lines have op- 
posite charge and collide, producing two additional vor- 
tex lines splintering away from the interaction area [see 
Fig- Btb)], but conserving the total vortex charge. The 
algebra of these vortex lines and the exciting possibility 
of controlling their direction will be discussed elsewhere. 

Stability of the planar vortex pairs depends critically 
on the coupling between vortices. We find no examples of 
stable out-of-phase co-rotating and the counter-rotating 
pairs, whereas the in-phase co-rotating vortices can be 
very robust. The in-phase configuration, due to the pos- 
itive coupling between nearest-neighbor lobes, is equiva- 
lent to a 7r-phase jump between the wavefunctions, due to 
the nature of the linear Bloch waves. Such phase "twist- 
ing" has a stabilizing effect as previously discovered in 
the framework of discrete nonlinear systems with attrac- 
tive nonlinearities ^3 • The actual region of stability for 
the co-rotating vortex states depends on the strength of 
the coupling (which depends on the values of \i and Vb), 
in agreement with the results for 2D vortices with an at- 
tractive nonlinearity . The state shown in Fig. 0Ja) 
is dynamically stable to t = 10 3 . 

Most remarkably, we find that vortex "stacks" of mul- 
tiple co-rotating vortices (somewhat analogous to the 
stacks of "pancakes" [l^). forming an extended discrete 
vortex line, can also be highly robust, in both the off-site 
[Fig. Eta)] and on-site configurations [Fig. ED 5 )]' These 
results demonstrate that a 3D optical lattice can provide 
transverse stabilization of vortex lines. This effect can 
be explored experimentally, since without the transverse 
periodicity, vortex lines exhibit strong 'snake-type' insta- 
bilities [12j • The possibility of long wavelength modula- 
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FIG. 5: (color online) Example of a stable off-site vortex 
stack at t = (a) and at t = fO 3 (b) (/i = 7.9). (c) Unstable 
on-site vortex stack at t = and t = f50 (fi = 7.5) (d). 

tional instability can not be ruled out however, especially 
at long time scales [lflj . Our numerical simulations per- 
formed with stacks up to 16 lattice sites long have found 
no trace of such an instability at t = 10 3 [Fig. EJa)]. 
When the instability is present, e.g. due to stronger 
inter-lobe coupling, it develops rapidly, leading to the 
spectacular "explosion" of the vortex line [see Fig. E{b)]. 

Apart from the generalizations of planar 2D gap vor- 
tices, 3D lattices can support vortices which have no 
2D analogs. The "truly" 3D vortex with the strongest 
coupling between lobes is the four-site state with simple 
phases (0, |, 1, §) (in multiples of it) shown in Fig.[|Jb). 
We term this configuration a "folded vortex" , which 
can be viewed as a combination of two dipoles (0, 1) 
and f ) that are not equivalent in inter-site coupling 
strengths. Similar coupling structure also occurs for a 
planar rhomboid vortex j2j . The simplest fully symmet- 
ric 3D vortex is that of the out-of-plane coupled dipole 
state of Fig. Etc). As can be seen by the particle flow 
arrows, the vortex circulation is highly nontrivial. 

Fully 3D lattices also offer us the possibility of com- 
bining different symmetry elements of 2D vortex states. 
The simplest example is the "vortex diamond" shown in 
Fig. Eld), which consists of two coupled planar vortices 
lying in the orthogonal symmetry planes. Unlike the in- 
coherently coupled out-of-plane vortices this is a co- 
herent fully 3D vortex characterized by a single chemical 
potential, and it can not be reduced to 2D equivalents. 

A remarkable feature of the uniquely 3D gap vortices 
is their dynamical stability up to at least t = 10 3 for 
a range of chemical potentials and lattice depths which 
is much broader than that of the bound states of two 
or more planar vortices. The experimental generation of 
these states is an open problem. However they appear 
to be an extremely robust means of introducing localized 
vorticity into the 3D lattice. 




FIG. 6: (color online) Novel 3D vortices at Vo — 6, [i — 7.9. 
(a) Folded vortex shown in (b) relative to the positions of 
the lattice minima, (c) Crossed-dipole vortex and (d) vortex 
diamond. Thick lines correspond to the primary vortex lines. 

In conclusion, we have predicted novel classes of spa- 
tially localized vortex structures in BECs loaded into 3D 
optical lattices. Such structures include planar vortices, 
vortex stacks, folded vortices and vortex diamonds. We 
have demonstrated numerically that many of these vor- 
tex structures are remarkably robust, and therefore they 
may be generated and observed in experiment. Our re- 
sults open up the possibility of studying more complex 
condensate structures such as localized vortex rings and 
knots in periodic potentials. 
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